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Abstract 

Background: Sphingosine kinase 1 (SKI) is a key regulator of the dynamic ceramide/sphingosine 1 -phosphate rheostat 
balance and important in the pathological cancer genesis, progression, and metastasis processes. Many studies have 
demonstrated SKI overexpressed in various cancers, but no meta-analysis has evaluated the relationship between SKI and 
various cancers. 

Methods:\Ne retrieved relevant articles from the PubMed, EBSCO, ISI, and OVID databases. A pooled odds ratio (OR) was 
used to assess the associations between SKI expression and cancer; hazard ratios (HR) were used for 5-year and overall 
survival. Review Manager 5.0 was used for the meta-analysis, and publication bias was evaluated with ST ATA 12.0 (Egger's 
test). 

Results: Thirty-four eligible studies (n = 4,673 patients) were identified. SKI positivity and high expression were significantly 
different between cancer, non-cancer, and benign tissues. SKI mRNA and protein expression levels were elevated in the 
cancer tissues, compared with the normal tissues. SKI positivity rates differed between various cancer types (lowest [27.3%] 
in estrogen receptor-positive breast cancer and highest [82.2%] in tongue squamous cell carcinoma). SKI positivity and high 
expression were associated with 5-year survival; the HR was 1.86 (95% confidence interval [CI], 1.18-2.94) for breast cancer, 
1.58 (1.08-2.31) for gastric cancer, and 2.68 (2.10-3.44) for other cancers; the total cancer HR was 2.21 (95% CI, 1.83-2.67; P 
< 0.00001). The overall survival HRs were 2.09 (95% CI, 1.35-3.22), 1.56 (1.08-2.25), and 2.62 (2.05-3.35) in breast, gastric, 
and other cancers, respectively. The total effect HR was 2.21 (95% CI, 1.83-2.66; P < 0.00001). 

Conclusions: SKI positivity and high expression were significantly associated with cancer and a shorter 5-year and overall 
survival. SKI positivity rates vary tremendously among the cancer types. It is necessary to further explore whether SKI might 
be a predictive biomarker of outcomes in cancer patients. 
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Introduction 

Sphingolipids are structural and functional components of 
biological membranes [1], and their metabolites, including 
sphinganine, ceramide, sphingosine, and sphingosine 1 -phosphate 
(SIP), have emerged as critical players in a number of fundamental 
biological processes. For example, these metabolites act as 
bioactive mediators in various cellular processes, including 
survival, proliferation, differentiation, and apoptosis [2,3]. More- 
over, sphingolipids are known to be involved in almost every type 
of disease [4] . They are reported to have regulation effect n cancer 
pathogenesis, progression, angiogenesis, proliferation, migration, 
inflammation, drug resistance, and cell death (apoptosis, necrosis, 
autophagy, and anoikis) [5,6]. Sphingosine kinases (SphKs), 
including sphingosine kinase type 1 (SKl/SphKl) and sphingosine 
kinase type 2 (SK2/SphK2), catalyze the phosphorylation of 
sphingosine to SIP and are crucial regulators of the balance 
among ceramides, sphingosine, and SIP [7]. Unlike SK2, SKI 
does not contain any transmembrane domains. In addition, their 



tissue distributions differ in human beings, as well as in animals [8— 
1 0] , and their opposing functions in sphingolipid metabolism have 
also been reported [11]. Specifically, SKI has become a cancer 
research hotspot in cancer and has recently been considered a 
bona fide oncogene [12]. 

Much evidence has shown that SKI can be detected in tumor 
tissues; notably, SKI was reported to be overexpressed in most 
studies [13-16]. Immunohistochemistry (IHC), polymerase chain 
reaction (PCR), and western blotting (WB) are usually used in such 
studies. The data showed that relative to paired non-cancer tissues, 
elevated SKI expression ranges from 2- to 8-fold at both the 
mRNA and protein levels in some cancers, including breast, lung, 
ovarian, stomach, and colon cancers. Meanwhile, some studies 
found that SKI overexpression might be related to cancer 
metastasis, reduced survival time, and poor prognosis [17,18]. 
Therefore, SKI was suggested as a novel biomarker of clinical 
prognosis in some cancers [19]- However, no systematic reviews or 
meta-analyses have discussed the role and clinical significance of 
SKI in cancer. Therefore, we performed a meta-analysis to 
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Figure 1. Flow chart of identified eligible studies. 
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investigate the relationship between SKI expression and various 
cancers and the possibility that SKI might be used as a cancer 
biomarker. 



the full texts were independently read and carefully checked to 
exclude ineligible studies. Disagreements were resolved via 
consensus. 



Materials and Methods 

Search strategy and study selection 

Electronic databases, including the English PubMed, EBSCO, 
ISI, OVID, ACS, and Cochrane Library databases and the 
Chinese VIP, Wan Fang, and China National Knowledge 
Infrastructure (CNKI) databases, and Google Scholar were 
searched for publications from inception to September 2013. 
The search keywords and terms were "sphingosine kinase type 1" 
OR "sphingosine kinase 1" OR "Sphkl" OR "Sphk-1" OR 
"SKI" OR "SK-1" OR "SK-I" AND "cancer" OR "tumor" OR 
"neoplasm" OR "carcinoma." Only Chinese- and English- 
language papers were included. All articles were imported into 
EndNote X6 software (Thomson Reuters Corporation, New York, 
NY, USA) to eliminate duplicates. 

All tides and abstracts were read by 2 independent reviewers 
(Yun ZHANG and Yan WANG) to select eligible studies. Next, 



Study selection criteria 

The inclusion criteria for primary studies were as follows: (1) 
human studies; (2) cancer types and detection method were clearly 
described; and (3) SKI expression was reported. 

The exclusion criteria for primary studies can be summarized as 
follows: (1) reviews, abstracts, and case reports; (2) animal or cell 
studies; (3) no description of detection methods; (4) a total sample 
size < 10; (5) a lack of SKI expression data or the IHC, PCR, and 
WB results were only reported in figures; (6) the hazard ratios (HR) 
and 95% confidence intervals (CI) of survival analyses were not 
reported or could not be calculated from other data; and all 
articles with the data from the same patient population, except for 
the most recent report. 

Data extraction 

A data collection form was designed before extracting data. The 
final included articles were assessed independendy by 2 reviewers 
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(Yun ZHANG and Yan WANG). The extracted data included the 
first author, year of publication, and country; cancer types, 
specimen types, and detection methods; the number of patients fn 
the experimental and control groups; the number of cancer, 
adjacent non-cancer, and benign samples and the number of SKI - 
positive (high expression or a score & ++) and negative patients 
(low expression or a score of - or +); the mRNA/protein 
expression levels in cancer and non-cancer tissues; the 5-year 
survival, overall survival (OS), disease-specific survival (DSS), and 
disease-free survival (DFS) outcomes and the follow-up duration. 
Among these, if the HR was not originally reported, we calculated 
or extracted this data according to the methods described by 
Tierney (2007) or Parmar (1998) [20,21]. 

Statistical methods 

The Review Manager software (version 5.02 for Windows; The 
Cochrane Collaboration, 2009) was used for the meta-analysis, 
and ST ATA 12.0 software (StataCorp, College Station, TX, USA) 
was used to analyze the publication bias (Egger's test). According 
to the Review Manager Handbook, for dichotomous data such as 
the number of patients with SKI positivity/high expression in 
cancer and non-cancer tissues, the odds ratios (OR) and 95% CI 
were combined to provide effective values. For continuous data 
such as SKI activity, the standardized mean difference (SMD) and 
95% CI were used to estimate effective values. For the association 
between SKI and OS, 5-year survival was described as a HR and 
95% CI. Statistical heterogeneity was examined with the 
Cochrane CHest (significant at P < 0.1) and I 2 value. When P 
< 0.1 or f < 50%, a fixed-effect model was used; otherwise, a 
random-effect model was used. If necessary, a sensitive analysis 
was also performed to evaluate the influences of individual studies 
on the final effect. When some studies were omitted or subgroup 
analyses were adopted and if no decreases in heterogeneity were 
observed, the qualitative systematic review method was used to 
describe the results. All P values were 2-sided, and P < 0.05 was 
considered significant. 

Results 

Search results and characteristics 

According to the inclusion criteria, 34 studies, including 1 1 
Chinese articles and 23 English articles, were eligible for the final 
meta-analysis (figure 1) [14-16,22-53]. The studies were published 
between 2005 and 2013. This analysis included 4,673 patients, 19 
types of cancer from 7 countries, and 5 types of detection methods. 
The follow-up duration ranged from 30 months to 12 years. Seven 
studies were about breast cancer; 6 studies, colon cancer; and 4 
studies, gastric cancer. The characteristics of the included studies 
are shown in Table 1. 
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SK1 expression in cancer and non-cancer tissues 

Initially, we assessed the association between SKI positivity/ 
high expression and cancer. Thirteen studies in which SKI 
expression was detected by IHC were included in this meta- 
analysis. Statistical heterogeneities were found in the subgroup 
analyses, except for that of adjacent non-cancer tissues versus 
benign tissues; therefore, a random-effect model was applied. The 
results demonstrated that significant differences in SKI positivity/ 
high expression rates were found between cancer and non-cancer 
tissues, cancer and adjacent non-cancer tissues, cancer and benign 
tissues, and adjacent non-cancer and benign tissues (all P < 
0.0001); the respective OR (95% CI) were 11.86 (6.37-22.08), 
6.66 (3.47-12.79), 11.94 (6.65-21.46), and 2.68 (1.70-4.20). SKI 
positivity/high expression is considered to associate with cancer 
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Figure 2. Forest plot of SKI positivity and high expression in cancer and on-cancer tissues. A subgroup analysis and random-effect 
model were used. The squares and horizontal lines correspond to the study-specific OR and 95% CI. The diamond represents the summary OR and 
95% CI. 
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(Figure 2). To prove the robust results regarding this association 
between SKI expression and cancer, a sensitivity analysis was 
conducted by omitting some of the obviously different studies until 



an acceptable heterogeneity level was reached (Table 2). In the 
absence of significant heterogeneity (P > 0.1, f < 50%) among 
the studies, a fixed model was adopted and the association was still 
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Table 2. Sensitivity analysis of high heterogeneity outcomes in meta-analysis. 
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Figure 3. SK1 mRNA and protein-level expression in cancers, compared with normal tissues, via RT-PCR (A) and western blot (B). 

SKI expression levels in all normal tissues are assumed to be 1. Histograms symbolize the relative levels of SK1 protein or mRNA expression in 
different cancers. 

doi:1 0.1 371 /journal.pone.0090362.g003 
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Table 3. SKI positive/high expression in variety cancers and 
non-cancer tissues. 



Cancer types 


Total 


SKI positive/high expression (%) 


ER+ breast cancer 


2202 


602(27.3) 


ER- breast cancer 


812 


477(58.7) 


Breast cancer(unknow) 


1287 


257(20.0) 


Breast cancer(total) 


4301 


1336(31.1) 


Salivary gland carcinoma 


159 


85(53.5) 


Esophageal cancer 


136 


110(80.9) 


Gastric Cancer 


412 


263(63.8) 


Colon cancer 


229 


170(74.2) 


Prostate cancer 


227 


163(71.8) 


Glioma 


243 


1 00(41 .2) 


Neuroblastoma 


142 


73(51.4) 


Thyroid Cancer 


42 


29(69.0) 


HNSCC 


18 


8(44.4) 


TSCC 


45 


37(82.2) 


Non-cancer tissues 


750 


255(34.0) 



Abbreviations SGC: salivary gland carcinoma; HNSCC: Head and neck 
squamous cell carcinoma; TSCC: Tongue squamous cell carcinoma; NA: not 
available. 

doi:1 0.1 371 /journal.pone.0090362.t003 



observed. Therefore, we think that our results were stable and 
credible. 

SK1 mRNA and protein expression level changes in 
cancers 

Seventeen studies reported SKI mRNA expression levels, and 
12 studies reported SKI protein expression levels. The sample 
sizes in these studies ranged from 4 to 95 cases. All outcomes 
indicated an increasing trend of SKI mRNA and protein 
expression levels in cancer tissues. When the SKI expression 
levels in normal tissues were assumed to be 1, the mRNA level 
increases ranged from 1.48-fold in uterine cancer to 7.86-fold in 
salivary gland carcinoma, and the protein level increases ranged 
from 1.45-fold in non-Hodgkin's lymphoma to 8.88-fold in gastric 
cancer (Figure 3). 



SK1 positivity/high expression rates in cancers 

Different cancer types present different rates of SK 1 positivity / 
high expression. We calculated the SKI positivity/high expression 
rates in various cancers, and these results are shown in Table 3. 
The SKI positivity/high expression rate in estrogen receptor- 
negative (ER— ) breast cancers (58.7%) was higher than that in 
ER-positive (ER+-) breast cancers (27.3%). The average positive 
rate in breast cancers was 31.1%. The lowest positive rate was 
found in breast cancer with unknown ER expression status, 
whereas the highest positive rate (82.2%) was found in tongue 
squamous cell carcinoma (TSCC; Table 3). Because of the small 
sample size of the TSCC study (n = 45), selection or other biases 
would unavoidably affect the result. 

SK1 enzyme activity in cancer tissues 

Four studies that measured SKI enzyme activity were 
evaluated. One was excluded because the SKI activity was 
measured in the blood [43]. A forest plot showed that the SMD 
(95% CI) were 3.27 (0.38-6.16), with P = 0.03 (Figure 4). One 
study by Leyre (2012) reported an obvious excess of SKI enzyme 
activity levels, compared with the other 2 studies [14]; the 
heterogeneity was so remarkable (P <0.00001, I 2 = 96%) that it 
could not be altered by omitting any single study from the 
sensitivity analysis; therefore, we cannot draw a firm conclusion 
regarding a relationship between SKI activity and cancer. 

Five-year and overall survival rates and SK1 expression in 
cancer tissues 

The HR values, log[HR], and SE[log(HR)] were either 
obtained directly from study data or extracted from survival 
curves, according to the methods described by Parmar or Tierney 
[20,21]. Generic inverse variance (IV) and subgroup analyses were 
adopted according to the types of cancer. As no significant 
heterogeneity was found, a fixed-effect model was used to perform 
the meta-analysis. The 5-year survival meta-analysis showed that 
in cancer patients, the HR of SKI positivity/high expression was 
1.86(95%CI, 1.18-2.94) for breast cancer (P = 0.008), 1.58(1.08- 
2.31) for gastric cancer (P = 0.02), and 2.68 (2.10-3.44) for other 
cancers (P < 0.0001). The total HR among all the cancers was 
2.21 (95% CI, 1.83-2.67; P < 0.00001). SKI positivity/high 
expression is thought to associate with a shorter 5-year survival 
duration in cancer patients (Figure 5). For overall survival, the HR 
was 2.09 (95% CI, 1.35-3.22) for breast cancer (P= 0.0009), 1.56 
(1.08-2.25) for gastric cancer (P = 0.02), and 2.62 (2.05-3.35) for 
other cancers (P < 0.00001). The total effect HR was 2.21 (95% 
CI, 1.83-2.66; P <0.00001). These data imply that SKI 



Cancer 



Non cancer 



Std. Mean Difference 



Study or Subgroup 



Mean SD Total Mean SD Total Weight IV. Random. 95% CI 



Std. Mean Difference 
IV. Random. 95% CI 



1.2.1 SKI activity in cancer(pmol/min/mg protein) 

PawetKnapp2010 8.11 0.9 23 3.11 0.38 

Leyre Brizuela 2012 259 55.2 13 88 75.4 

Bernard Malavaud 2010 11.1 8.4 28 5.9 3.2 
Subtotal (95% CI) 64 
Heterogeneity: Tau 2 = 6.19; Chi 2 = 48.28, df = 2 (P < 0.00001 ); I 2 
Test for overall effect: Z = 2.21 (P = 0.03) 



Total (95% CI) 

Heterogeneity: Tau 2 = 
Test for overall effect: 



64 



31.6% 
33.6% 
34.8% 
100.0% 



6.19; Chi 2 = 48.28, df = 
Z = 2.21 (P = 0.03) 



2 (P < 0.00001); r 



= 96% 



59 100.0% 

= 96% 



6.80 [5.13, 8.47] 
2.51 [1.44, 3.57] 
0.81 [0.26, 1.35] 
3.27 [0.38, 6.16] 



3.27 [0.38, 6.16] 



-10 -5 0 5 10 

Favours non-cancer Favours cancer 



Figure 4. Forest plot of the association between SKI enzyme activity and cancer tissues. A random-effect model was used. The pooled 
SMD is symbolized by a solid diamond at the bottom of the forest plot, the width of which represents the 95% CI. 
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Study or Subgroup 



loqfHazard Ratio] SE Weight 



Hazard Ratio 
IV. Fixed. 95% CI 



Hazard Ratio 
IV. Fixed. 95% CI 



1.3.1 Breast cancer 

Carol Watson 2010 0.683 0.35 

Eugen Ruckhaberle 2013 0.255 0.98 

Jan Ohotski 2013 0.61 0.33 
Subtotal (95% CI) 

Heterogeneity: Chi 2 = 0.17, df = 2 (P = 0.92); I 2 = 0% 
Test for overall effect: Z = 2.67 (P = 0.008) 



7.6% 
1.0% 
8.5% 
17.1% 



1.98 [1.00, 3.93] 
1.29 [0.19, 8.81] 
1.84 [0.96, 3.51] 
1.86 [1.18, 2.94] 



1.3.2 Gastric cancer 

Wen Li 2009 0.31 0.24 16.1% 1.36 [0.85, 2.18] 

Zhuge Yonghu 2011GC l-ll 0.723 0.37 6.8% 2.06(1.00,4.26] 

Zhuge Yonghu 2011GCIII-VI 0.751 0.69 1.9% 2.12(0.55,8.19] 

Subtotal (95% CI) 24.8% 1.58 [1.08, 2.31] 
Heterogeneity: Chi 2 = 1 .07, df = 2 (P = 0.58); I 2 = 0% 
Test for overall effect: Z = 2.37 (P = 0.02) 



1.3.3 Other cancer 

Guan Hongyu 2009 
Jian Pan 2011 
Liu Guanglin 2010 
Subtotal (95% CI) 

Heterogeneity: Chi 2 = 0.17, df = 2 (P = 0.92); I 
Test for overall effect: Z = 7.82 (P < 0.00001 ) 



1.082 
0.956 
1.061 

: 0% 



0.36 
0.15 
0.31 



7.2% 
41.3% 

9.7% 
58.1% 



2.95 [1.46, 5.98] 
2.60 [1.94, 3.49] 
2.89 [1.57, 5.30] 
2.69 [2.10, 3.44] 



Total (95% CI) 

Heterogeneity: Chi 2 = 7.37, df = 8 (P = 0.50); I 2 = 0% 
Test for overall effect: Z = 8.24 (P < 0.00001 ) 

Test for subgroup differences: Chi 2 = 5.96, df = 2 (P = 0.05), I 2 = 66.4% 



100.0% 2.21 [1.83, 2.67] 



0.1 0.2 0.5 1 2 5 10 

Favours ski low/negative Favours ski high/positive 



Figure 5. Forest plot of association between SKI expression and 5-year survival. A subgroup analysis for different cancers and fixed-effect 
model were used. The squares and horizontal lines correspond to the study-specific HR and 95% CI. The areas of the squares reflect the weights 
(inverse of the variance). The diamond represents the summary HR and 95% CI. 
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positivity/high expression is related to overall survival in various 
cancers (Figure 6). 

Publication bias 

Egger's test was performed to assess the publication bias in the 
literature. In a subgroup analysis of SKI expression in cancer, 
adjacent non-cancer, and benign tissues, publication bias was 
found for cancer versus adjacent non-cancer (P = 0.009) and 
cancer versus benign tissues (P = 0.04); the funnel plot was not 
symmetrical. No publication bias was found for cancer versus non- 
cancer (P= 0.059) and adjacent non-cancer versus benign tissues 
(P = 0.176; Figure 7). This might be a limitation of our analysis 
because studies with null findings, especially those with small 
sample sizes, were less likely to be published. In addition, different 
test methods were used in these studies. No publication bias was 
found with regard to 5-year (P = 0.754) and overall survival 
(P = 0.175; Figures 8 and 9). 

Discussion 

The lipid kinase SKI is expressed extensively in human beings, 
mice, yeast, and plants and plays an important role in regulating 
the dynamic balance of the ceramide/SIP rheostat. SKI 
overexpression or increased activity leads to SIP accumulation 
that in turn contributes to prosurvival and anti-apoptotic 
mechanisms, as well as cancer genesis, metastasis, and drug- 
resistance. In contrast, ceramide accumulation leads to pro- 
apoptotic, autophagic, and anti-survival effects. High SKI 



expression was detected in many cancers, including those of the 
lung, breast, ovary, stomach, and kidney. 

IHC of formalin-fixed, paraffin-embedded samples is among the 
most commonly used methods to detect SKI expression. We 
evaluated the association of SKI and different cancers and non- 
cancer tissues according to the IHC outcomes. In the subgroup 
and total effect analyses, SKI positivity/high expression was 
significandy increased in cancer tissues, compared with non- 
cancer (adjacent and/ or benign) tissues. The OR between cancer 
and adjacent tissues (2.68) was less than that between cancer and 
benign tissues (11.94). This finding demonstrates that the SKI 
expression level gradually increased in benign, adjacent, and 
cancer tissues. This might indicate that SKI expression is 
associated with degrees of pathological differentiation or hetero- 
geneity in cancer. However, different IHC judgment criteria 
existed among these studies. Positive and negative were used to 
judge some, while others were graded according to signs (— , +, ++, 
and +++). In our meta-analysis, we used negative/low (including 
— and +) and positive/high (including ++ and +++) to assess SKI 
expression in different tissues. That might be partly responsible for 
the heterogeneity. 

When the mRNA and protein expression levels were detected 
via reverse transcriptase PCR or western blotting, the results were 
consistent. We found that the SKI positive/high expression rates 
varied tremendously between different types of cancer (ranging 
from 20-82.2%). Differences in SKI tissue distribution might be 
the main reason for this variation. Because the sample sizes and 
detection methods were different, the high heterogeneity prevent- 
ed us from merging the studies to perform a meta-analysis. 



PLOS ONE | www.plosone.org 



8 



February 2014 | Volume 9 | Issue 2 | e90362 



Sphingosine Kinase 1 and Cancer 



Study or Subgroup 



lo qfHazard Ratio] SE Weight 



Hazard Ratio 
IV. Fixed. 95% CI 



Hazard Ratio 
IV. Fixed. 95% CI 



1.4.1 Breast cancer 

Eugen Ruckhaberle 2013 0.215 0.83 

Jan Ohotski 2013 0.775 0.23 
Subtotal (95% CI) 

Heterogeneity: Chi 2 = 0.42, df = 1 (P = 0.52); I 2 = 0% 

Test for overall effect: Z = 3.32 (P = 0.0009) 

1.4.2 Gastric cancer 

Wen Li 2009 0.3 0.23 

Zhuge Yonghu 2011GC l-ll 0.723 0.37 

Zhuge Yonghu 201 1GCIII-VI 0.751 0.69 
Subtotal (95% CI) 

Heterogeneity: Chi 2 = 1.16, df = 2 (P = 0.56); I 2 = 0% 

Test for overall effect: Z = 2.35 (P = 0.02) 



1.4.3 Other cancer 

Brianna N. Harris2013 -0.236 0.89 

Guan Hongyu 2009 1.082 0.36 

Jian Pan 2011 0.956 0.15 

Liu Guanglin 2010 1.061 0.31 
Subtotal (95% CI) 

Heterogeneity: Chi 2 = 2.03, df = 3 (P = 0.57); I 2 = 0% 
Test for overall effect: Z = 7.71 (P < 0.00001) 



Total (95% CI) 

Heterogeneity: Chi 2 = 9.03, df = 8 (P = 0.34); I 2 
Test for overall effect: Z = 8.40 (P < 0.00001 ) 
Test for subgroup differences: Chi 2 = 5.43, df = 



11% 



1 .3% 
16.8% 
18.1% 



16.8% 
6.5% 
1.9% 

25.2% 



1.1% 
6.9% 
39.5% 
9.3% 



1.24 [0.24, 6.31] 
2.17 [1.38, 3.41] 
2.09 [1.35, 3.22] 



1.35 [0.86, 2.12] 
2.06 [1.00,4.26] 
2.12 [0.55, 8.19] 
1.56 [1.08, 2.25] 



0.79 [0.14, 4.52] 
2.95 [1 .46, 5.98] 
2.60 [1.94, 3.49] 
2.89 [1.57, 5.30] 



56.7% 2.62 [2.05, 3.35] 



100.0% 2.21 [1.83, 2.66] 



2(P = 0.07), I 2 = 63.1% 
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Figure 6. Forest plot of association between SKI expression and overall survival. A subgroup analysis for different cancers and fixed-effect 
model were used. The squares and horizontal lines correspond to the study-specific HR and 95% CI. The areas of the squares reflect the weights 
(inverse of the variance). The diamond represents the summary HR and 95% CI. 
doi:1 0.1 371 /journal.pone.0090362.g006 




Figure 7. Funnel plot for publication bias in an association analysis of SK1 expression among cancer, adjacent non-cancer, and 
benign tissues. Begg's linear regression test was used, and P = 0.059, 0.009, 0.040, and 0.176 in a subgroup analysis of SKI positivity and high 
expression in cancer and non-cancer tissues, respectively (Figure 2). 
doi:1 0.1 371 /journal.pone.0090362.g007 
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Begg's funnel plot with pseudo 95% confidence limits 
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Figure 8. Funnel plot for publication bias regarding 5-year survival. Begg's linear regression test was used; P = 
doi:1 0.1 371 /journal.pone.0090362.g008 



0.754. 



Moreover, only 3 studies [14,28,30] reported SKI enzyme 
activity, and the interstudy activity level gaps were very wide 
(Figure 3). In all the 3 studies, the radioassay method was used to 
detect SKI activity. This gap might have been introduced by the 
operator or testing instrument. It might also have been caused by 
differences in cancer types. To obtain more accurate results, more 
studies and the evaluation of more samples of the same cancer type 
with identical detection methods are needed. 

Regarding the survival analysis, the follow-up duration in these 
studies ranged from 30 months to 12 years. We used 5-year and 
overall survival to evaluate the correlation between SKI expres- 
sion and survival in cancer patients. Studies in which the follow-up 



time was S5 years were included in the 5-year survival analysis. 
Because the overall survival follow-up times differed, we could not 
acquire concrete data from some of the articles and instead 
extracted data from survival curves. Therefore, the outcomes 
might not be very accurate. However, we do not believe that this 
affected the total outcome. 

In our study, SKI expression was detected in 19 common types 
of cancer. Some of the more common cancers were the focus of 
more studies, including breast cancer, which was researched in 7 
studies [23,25,31,32,44,51,53], and colon cancer, which was 
researched in 5 studies[15, 27, 45,49,52]. Therefore, different 
cancers were represented unequally with respect to sample 



Begg's funnel plot with pseudo 95% confidence limits 
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Figure 9. Funnel plot for publication bias regarding overall 

doi:1 0.1 371 /journal.pone.0090362.g009 



survival. Begg's linear regression test was used; P = 0.175. 
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numbers. Meanwhile, the cancer tissue sample and control sample 
selection methods differed among different studies. In some 
studies, the cancer and control tissue samples were collected from 
different areas in the same patient, whereas in other studies, the 
samples were collected from different patients. Moreover, in these 
studies, different SKI detection methods were used and cancer 
patients of different ethnicities were included. All of these could be 
sources of heterogeneity. Despite the heterogeneity that was 
observed in some subgroup analyses, the total outcomes were not 
affected by the sensitivity analysis. 

There were some other limitations and potential biases in our 
meta-analysis. First, as few (<3) studies regarding disease-specific 
survival (DSS) and disease-free survival (DFS) were included, we 
could not perform a meta-analysis of those variables, which are 
used to evaluate prognosis in cancer patients. Second, this meta- 
analysis did not discuss the associations between SKI and 
ethnicity, tumor size, TNM stage, histological classification, 
pathological differentiation, lymphatic invasion, and metastasis, 
all of which are important clinical characteristics in cancer. In 
addition, only 1 or 2 studies were identified for some types of 
cancer, and as those sample sizes are insufficient, the associated 
biases might be remarkable. 

In the individual papers, the SKI expression were found in one 
or two kinds of cancers in each study. In our meta-analysis, we 
found that there were relationships between SKI expression and 
various cancers, which may imply the possibility of using the SKI 
level as the clinical biomarker of cancers. 
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